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Abstract 

> 

0^ ' The latest structure function results, as presented at the ICHEP98 conference, are 

. reviewed. A brief introduction to the formalism and the status of global analyses of parton 

\ distributions is given. The review focuses on three experimental areas: fixed-target results 

and their constraints on the parton densities at high x; spin structure and spin parton 
■ densities as well as the status of the associated sum rules; HERA results on the dynamics 

\ of F2 at low (x, Q^), charm and Fl as well as the measurement and interpretation of the 

Q \ high-Q^ neutral and charged current cross-sections. 

i 
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^ ■ Slides are available from 



O . Plenary talk presented at the XXIX ICHEP98 Conference, 

Vancouver, July 1998. 



http://www-zeus. desy.de/conf erences98/#ichep9i 



1 Introduction - Formalism and Road Maps 



The differential cross-section l{k)N{p) — » l{k')X{p') for a lepton (e, /i) with four-momentum k 
scattering off a nucleon with four-momentum p can be expressed as 

(Pa{l^N) _ 2Txa^ 
dxdQ^ ~ xQ^ 

■ [F+F2(x, Q2) ^ Y_xF^{x, Q2) _ y-^FL{x, Q^)] 

where is the four-momentum transfer squared, x = Q'^/2p-q is the Bjorken scahng variable, 
y = p ■ q/p ■ k is the inelasticity variable and Y± = 1 ± (1 — y)^. The contribution from F2 
dominates the cross- sect ion. The contribution from Fl is a QCD correction which is important 
only at large y and that from xF^ is negligible for « M|. To investigate sensitivity to 
Fl at large y or XF3 at large Q^, the reduced cross-section a = ^^y^ dfdQ'^ adopted. In 
the Quark Parton Model (or in the DIS scheme of NLO QCD) and for Q'^ « M| F2/X is the 
charge-weighted sum of the quark densities 

F2(a;,g2)=a;^e2-S(x,Q2) 

i 

where E(x, Q^) = '}li[(li{x,Q'^) + Q^)] is the singlet summed quark and anti-quark dis- 
tributions. Similarly, the charged current cross-section e^{e~)N v{i')X at HERA can be 
expressed as 

dxdQ'^ 4nx\Q^ + M^) 

■[Y+F2^^{x,Q')tY.xF^''{x,Q% 

For fixed-target z/(z/)A^ /i+(/i^)X experiments, « M^, and in the QPM 

xF^''{x,Q^) = xY,e^-qNs{x,Q') 

i 

where qNsi^, Q"^) = J^ilQii^: Q"^) ~ <li{x, Q"^)] is the non-singlet difference of these distributions. 

In Fig. 0, the kinematic plane covered by the F2{x,Q'^) measurements is shown, including the 
new preliminary datasets from HI and ZEUS which are seen to extend to: low y (?/hera ~ 0.005) 
providing overlap with the fixed-target experiments; very low x (x^lO~^) at low exploring 
the transition region from soft to hard physics; high y {y ^ 0.82) giving sensitivity to F^, 
high X — >■ 0.65 probing sensitivity to electroweak effects in F2 and xF^ as well as constraining 
the valence quarks at large Q^. The fixed-target experiments NMC, BCDMS, E665 and SLAG 
experiments have provided final measurements at higher x and lower Q^. New information was 
presented at ICHEP98 from CCFR, E866 and the Tevatron which also constrain the medium- 
high X partons. 

Theoretically, the directions in [x, Q"^) can be mapped out according to the dominant dynamical 
effects, as illustrated in Fig. 0. Given a phenomenological input as a function of x, the parton 
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distributions are evolved to different physical scales (Q^) via the DGLAP evolution equations. 
The alternative BFKL approach is to attempt to calculate the x dependence directly from 
pQCD, where the running of the effective coupling constant is neglected to leading order. BFKL 
predicts an x"'^ dependence of F2 at small x. The BFKL equation has recently been calculated 
to NLO. These corrections are numerically very large in the experimentally accessible low 
(x, Q^) range, resulting in A^o ~ 0.5 being reduced to Xnlo ~ 0.1. The understanding of the 
NLO corrections to the BFKL equation is therefore an active area of study. In the DLLA 
('double leading log approximation'), non-leading In(Q^) ln(l/x) terms can also be evaluated, 
but a method which reliably maps the complete region of (x, Q^) in terms of pQCD is still 
not known. The region of high parton density may be reached at very low x, where these 
approaches are not strictly valid. The expectation from the GLR equation is that the region 
where the partons overlap is accessible at slightly higher x for decreasing Q^. However, 
should also be sufficiently large that higher twist and non-perturbative effects parameterised in 
terms of Regge exchanges can be neglected. 

In the DGLAP approach, the non-singlet contribution evolves as 

dqNs{x,Q^) _ a,(Q^) .yg ^ 
di ^qNs{x,Q) 

where t = ln((5^/A^) and the Pi/s represent the NLO DGLAP splitting probabilities for radi- 
ating a parton with momentum fraction x from a parton with higher momentum x'. Quantities 
such as xFs provide a measure of as(Q^) which is insensitive to the a priori unknown gluon 
distribution. Similarly, the singlet quark and gluon densities are coupled via 

/S(x',Q2)\ 
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\g{x',Q 

and quantities such as F2 provide input for S(x,(5^) as well as coupled knowledge of ^^(Q^) 
and the gluon, g{x,Q^). 

At the starting scale, Ql, the light valence quarks (g„) and the sea of quark and anti-quarks 
{q) as well as the gluon [g) are attributed a given functional form. For example in the MRST 
parameterisations 

/,(x, Ql) = Aix'^il - xr^il + e,v^ + ^,x) 

where some parameters are set to or fixed by sum rules and differences of u and d quarks may 
additionally be constrained. The heavy quark (sea) contributions are calculated explicitly at 
NLO and their uncertainty is typically determined by the range allowed for the effective mass 
of the heavy quark. The measured structure functions are then described by the convolutions 
of the light quark densities with the appropriate NLO matrix elements. 

The outlined procedure defines the structure of a nucleon in terms of its constituent quarks 
and gluons. However, guided by the new datasets and analyses which were presented to this 
conference, it should be noted that the following assumptions are made: 

• a;s(Q^) In(Q^) (DGLAP) terms are large compared to a<((Q^) ln(l/x) (BFKL) terms in 
the perturbative splitting functions; Q 



The study of inclusive quantities such as F2 at smah x are presently unable to distinguish these BFKL 
terms. The status of forward jet production searches which enhance sensitivity to these effects is reviewed by 
J. Huston. 01 
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Figure 1: Measured regions of F2 in the (x, Q^) kinematic plane. The nominal acceptance region 
of the HERA measurements corresponds to ?/hera > 0.005. The fixed-target experimental data 
occupies the region of high x at low Q^. 




Figure 2: Schematic representation of parton densities and the theoretical evolution directions 
in the (x, Q^) kinematic plane. 
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Figure 3: MRST parton distributions xf{x,Q^) at = 20 GeV^ 

• higher-twist (HT) contributions (suppressed by factors of 1/Q^) are neghgibly small; Q 

• Nuclear binding effects are treated as small corrections or ignored in analyses of deuteron 
data. 

In Table |l], the experimental datasets considered in the MRST analysis are listed along with the 
underlying physics process and the parton behaviour which is being probed. The experiments 
denoted by f correspond to final measurements reported at ICHEP98. Those denoted by 
correspond to new preliminary measurements reported at this conference. The latest global fits 
of MRST 0, GRV98 [| and CTEQ4 are used here to compare with the data. 

The output MRST parton distributions evolved to = 20 GeV^ are shown in Fig. 0. At this 
Q^, the expected fractions of the total momentum carried by the valence quarks is 25% and 

2) Progress in determining the size of these higher twist and other hadronisation power corrections via infra-red 
renormalons is reviewed by Y. Dokshitzer. ||^ 

•'^A discussion of diffractive final states in D.I.S including nuclear effects is included in the review by M. 
Erdmann. |^ 
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10%(4)- 6%(2u), 8%{2d), 5%(2s) and 2%(2c) is carried by the qq symmeteric sea and U%{g) 
is carried by the gluons. 

Incorporating new datasets and theoretical understanding improves the precision with which 
individual parton densities are known. In particular, uncertainties in d/u at high x, d/u at 
intermediate x, s at all x values, g at low x and c at low x, are discussed in the context of the 
presentations made at ICHEP98. 



2 Fixed- Target Results 
2.1 Determination of a. 



Gross-Llewellyn Smith (GLS) Sum Rule: The GLS sum rule expresses the fact that there are 
three valence quarks in the nucleon, subject to QCD corrections 



/'xF3(x,g^)- = 3(l-^-a2(^) 



,ac,n, AHT 



This is a fundamental prediction of QCD which, as a non-singlet quantity, is independent of the 
gluon distribution. The sum rule has been calculated to 0{a^) and estimated to 0{a^). [1^, |ll 



The CCFR collaboration have now published their final results, [0 incorporating earlier neu- 
trino measurements (WA25, WA59, SKAT, FNAL-E180 and BEBC-Gargamelle) to determine 
the GLS sum rule. In Fig. |^, the sum rule is proportional to the area under the data in four 
regions of Q^. The analysis of the world data enables as{3 GeV^) to be determined at NNLO 
accuracy and evolved to M| as 

The largest contribution to the systematic error is the uncertainty on the ratio of the total neu- 
trino and anti-neutrino cross-sections, (Tp/a,^ = 0.499 ±0.007 which determines the overall nor- 
malisation of xF^. This uncertainty will be improved with the NuTeV tagged u-u beamline. The 



nuclear target corrections to the GLS sum rule are predicted to be small. [|T^ The largest theory 
uncertainty is that associated with the higher twist contribution, AHT = 0.15±0.15 GeV^. 
Here, a renormalon approach, where chains of vacuum polarisation bubbles on a gluon propaga- 
tor lead to a prediction of the higher twist contribution in the perturbative expansion, predicts 
a small correction, AHT < 0.02 GeV^. Other models predict significantly larger HT corrections 
and the uncertainty encompasses this range. The renormalon approach is now rather successful 
in describing a range of power corrections to hadronic final state variables and leads to a 
central value of as{M'^) = 0.118. The most recent analysis of the xF^ data determines the 
same central value with a similar theoretical uncertainty. 



Scaling Violations at Large x: The most accurate method to determine as(M|) from the 
CCFR data remains the measurement of the scaling violations of the structure functions, 
F2 and xFs, using a NLO QCD fit. |jl5[ The scaling violation slopes (i(logF2)/(ilogQ^ and 
d{\og xF3)/d log Q"^ are shown in Fig. |^. From F2, the high-x gluon is also constrained to be 



xg{x, Ql = 5 GeV^) = (2.22 ± 0.34) x (1 - x 



,4.65±0.68 
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Figure 4: CCFR analysis of xF^ versus logio{x) for four ranges of Q^. The curve indicates 
a power law (Ax^) fit applied for x < 0.1 used to determine the integral in the unmeasured 
shaded region. 
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in agreement with global analyses using prompt photon data. From the combined results on 
F2 and xF^ 

asiMl) = 0.119 ± 0.002(exp.) ± Omi{HT) ± 0.004(sca/e) 

which represents one of the most precise determinations of this quantity. This improved mea- 
surement is higher than the earlier CCFR value due to the use of a new energy calibration. 
It is also higher than the SLAC/BCDMS analysis value of = 0.113 ± 0.003(exp.) ± 

0.004{theory). [jl6[ Here it is noted that there is a small but statistically significant discrepancy 
between the SLAC and BCDMS data (see the x > 0.45 points in Fig. which can be 
resolved if the correlated systematics of the BCDMS data at low y are taken into account. The 
quoted central value of the combined analysis, however, does not take these systematics into 



account. |T8[ There is therefore little evidence for any discrepancy with respect to the world 
average value af^^(M|) = 0.119 ± 0.002. 



Improved statistics will enable the NuTeV collaboration to determine decoupled from the 
gluon using the xF^ data alone. Further progress in the F2 analysis requires the calculation of 
NNLO terms in order to reduce the renormalisation and factorisation scale uncertainties. In 
this way the overall uncertainty on q;s(M|) can be reduced to ±0.002. In this area, there is 
still an outstanding 20% discrepancy between the CCFR F2 and NMC F2 data in the region of 
X ~ 10~^ which is discussed in the context of the strange quark sea and the ZEUS preliminary 
data later. However, this effect is negligible in the determination of a^, which depends mainly 
upon the high-x data as seen in Fig. ||. 



2.2 Nucleon Structure 

d/u at large x: Valence parton distributions at high x received attention at this conference 
provoked by new data and corrections of the d/u ratio from an analysis by U. Yang and A. 
Bodek. Extractions of from NMC and SLAC deuteron data have previously accounted 



for Fermi motion but not nuclear binding effects. A physically appealing model by Frankfurt 



and Strikman pO[ assumes that binding effects in the deuteron and heavier nuclear targets 
scale with the nuclear density. This nuclear binding correction is about 4% at x = 0.7 for 
fixed-target deuteron experiments, which is parameterised as an additional term 6{d/u) = 
(0.10±0.01)x(l+a:) added to the MRS(R2) PDF. This is sufficient to increase the d distribution 
significantly at high = 10^ GeV^ by about 40% at a; = 0.5, due to DGLAP evolution of the 
partons. The modification gives an improved fit to the NMC deuteron data. It is compelling 
in that this simple modification now improves the description of high-x CDHSW up/ up data 
(not shown) as well as the new data on the W asymmetry from CDF (probing intermediate 
X ~ Mw/y/s ~ 0.05 values) and the charged current cross-sections from ZEUS (see later), data 
which are free from nuclear binding effects. 

The world and F2 fixed-target data at high x is plotted in Fig. Here the NLO pQCD cal- 
culation incorporates target mass (TM) effects determined using Georgi-Politzer scaling where 



the scaling variable x is replaced by ^ = 2x/(l + yl + 4:M'^x'^/Q'^) as well as the correction 
for nuclear binding effects for the deuteron data. The description of the low data is signif- 
icantly improved {x^/DoF = 1577/1045) if higher twist (HT) corrections are incorporated as 
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indicated by the fit denoted NLO(pQCD+TM+HT). |jT9| A good description is also obtained 
for tlie very liigli x SLAG DIS data (not sliown) up to x = 0.9, altliougli tliere is some residual 
dependence in the resonance region (0.9 < x < 1). The renormalon approach therefore 



successfully predicts the x dependence of the higher twist terms although the normalisation 
constant is half of the previously estimated value. This renormalon approach also successfully 
describes the dependence of R = <7l/(Tt (not shown). In a NNLO fit this higher twist 
contribution tends to zero which suggests that the non-singlet higher twist contribution (rel- 
evant for xF^) is also small, as discussed in relation to the GLS sum rule. In addition, new 
measurements of R, covering the range of 0.015 < x < 0.55 from CCFR were presented which 
are in agreement with those from EMC, BCDMS and SLAC and extend the kinematic range 
up to g2 ^ 100 GeV\ 

As can be seen in Fig. ^, the u quarks in the proton carry more momentum on average than 
the d quarks. (The d quarks carry only slightly more momentum than the u quarks at high 
X.) and W~ production in pp collisions is primarily due to the annihilation processes 

u{p)d{p) — > and u{p)d{p) — >• W~ . On average the W^'s are therefore produced relatively 
more forward than the W~^^. CDF measure the lepton asymmetry from the W decays. 



da^/dyi — da jdyi 
da+ /dyi + da~ /dyi ' 



where da^ /dyi are the cross-sections for decay leptons as a function of lepton rapidity 
relative to the proton beam direction. The W decay to leptons tends to reduce the asymmetry 
at forward \yi\ but in a well-defined manner (assuming the W decays proceed via a SM V-A 
interaction). In Fig. |^, the final results from runs la and lb are shown compared to various 
PDF's using the DYRAD NLO calculations. The best agreement with the data is obtained 
using the modified MRS(R2) or the MRST PDF (not shown), which incorporates this data as 
input, indicating that the ratio of (i/w is larger than previously thought. 

d/u at intermediate x: The structure function measurements constrain the sum of d and u but 
not the difference. The classic method to constrain this difference is via the Gottfried sum rule 
where 

Sg = Fi- F^— = - {d,- u,)dx - - / {d-u)dx. 
Jo X 6 Jo 6 Jo 

The NMC result for 0.004 < x < 0.8 extrapolated over all x gives Sq = 0.235 ± 0.026, 
significantly below 1/3. [^] It should be noted however that the deuterium data used to derive 
F2 are not corrected for the effects noted above. This correction is estimated to be —0.013 i.e. 
50% of the uncertainty. |^ 



At this conference, the E866 data were reported, which provide a significant constraint 



on the ratio of d/u via the measurement of Drell-Yan dimuons with mass M^+^- > 4.5 GeV 
from 800 GeV protons on proton and deuterium targets as shown in Fig. |^. Here, the ratio of 
a^'^/2a^^ is measured as a function of the reconstructed momentum fraction of the target quark 
in the range 0.036 < X2 < 0.312. The MRST calculations incorporate a direct parameterisation 
of d — u in order to fit the data. Earlier PDF's such as MRS(R2) clearly overestimate this 
asymmetry but a, d = u sea is still ruled out. 
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1 5 10 50 100 1 5 10 50 100 



Q2 [Proton] Q2 [Deutron] 

Figure 6: Comparison of high-x and F2 data versus fitted with and without higher 
twist (HT) corrections. The error bars correspond to statistical and systematic errors added in 
quadrature. 
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Figure 7: CDF lepton charge asymmetry from W decays as a function of lepton rapidity. The 
error bars correspond to statistical and systematic errors added in quadrature. 



The E866 data require an overall decrease of the sea compared to the MRS(R2) parameter- 
isation. This is seen in the lower part of Fig. ^ where the E866 data are plotted as the sea 
contribution to -Ffl^^) ~ -^2* (2^)? which contributes negatively to So, evolved to the NMC 
values. The NMC data are plotted in the upper plot, compared to the valence as well as the 
summed contribution for the MRST and MRS(R2) partons. Although the valence and sea 
distributions differ, these distributions give similar results for the total -Ff (^) ~ F2{x). How- 
ever, it is clear that the agreement of the NMC data with either parameterisation is poor at 
intermediate x, which would require additional changes to the valence quark distributions. The 
MRST NLO fits determine So = 0.266, above the NMC measured value discussed earlier. This 
discrepancy is induced in part by the inclusion of the E866 data in the fit. 

A further constraint is provided via the HERMES semi-inclusive measurement of the ratio of 
the differences between charged pion production for (unpolarised) proton and neutron (^He) 
targets 



r(x, z) 



Nf{x,z)-N^'{x,z) 
N;^{x,z)-Nf{x,zy 



also reported at this conference. p8[ The ratio is observed to scale as a function of the frag- 



mentation scaling variable z and the data are used to constrain {d — u)/{u — d) at leading order 
as a function of x, as shown in Fig. p!0|a. Using the PDF's to constrain the valence quarks, the 
d — u distribution obtained from the E866 and (lower Q"^) HERMES data are shown to be in 



agreement in Fig. 10 b 



Generically, the Pauli exclusion principle, suppresses production of sea uu quarks at large x due 
to the additional quark. Hence the excess of d over -u is a direct consequence of the excess of 
over d^. The currently favoured dynamical approach is to include non-perturbative effects 
from virtual mesons directly from |n7r"'" > Fock states in a virtual pion model or to derive the 
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Figure 8: E866 ratio of 0-^^/2(7^'^ versus X2, the momentum fraction of the target quark, com- 
pared to the NLO calculations discussed in the text. 
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Figure 9: NMC Ff' — versus x compared to LO predictions discussed in the text. For each 
prediction, the top (bottom) curve is the valence (sea) contribution and the middle curve is the 
sum. The E866 results for the sea quark contribution to F|' — F2 are indicated as the negative 
contribution. 
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Figure 10: a) [d — u){u — d) extracted directly from HERMES data (•) compared to various 
parameterisations of the proton, b) HERMES data converted to — u in order to compare with 
E866 (o) data. 



sea from the valence quarks coupling to Goldstone bosons (e.g. vr"*" from u —* dir^) in a chiral 



model. ||2^ In conclusion, the latest data indicate that the ratio of d/u is less than previously 



thought, a decrease which is correlated with the increase of d/u. 

s quark contribution: The strange contribution to the quark sea is determined using opposite- 
sign dimuon events in uFe scattering from CCFR where one muon comes from the lepton 
vertex and the other results from the semi-leptonic decay of a charmed hadron. The ratio of s 
to u and d quarks integrated over x in an NLO analysis was determined to be 0.477lo;o53. 



This large uncertainty requires, for example, that NuTeV tag z/ and z/ and take the difference 
of these cross-sections in order to reduce the uncertainty on sin?6w, also reviewed at this 



conference. A second potential method to constrain the strange quark contribution is by 
comparing Fg^^ and Flf^ . The analysis of this data was discussed by C. Boros. To leading 
order the muon (NMC) and neutrino (CCFR) structure functions are related by the "5/18ths 
rule" 



piN _ _ 3 a^g(a^) + xs{x) ^pj^j^ 
^ ~ 18^ bxq{x) +xq{xy ^ 
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where the [ss symmetric) strange sea enters as a correction. However, the CCFR data cor- 
rected using the dimuon result to constrain the strange sea hes significantly (~ 20%) above 
the NMC data for all values of at x < 0.1. Here nuclear shadowing corrections, estimated 
independently for uFe (CCFR) and /iD (NMC) data, indicate that part of the discrepancy 



between the two experiments at low x may be accounted for in this way. |^ However, there 
remains a significant discrepancy: the strange quark sea and the relation between and ^2'^ 
data for x < 0.1 is therefore an active area of theoretical and experimental investigation. 



3 Spin Structure Functions 

Major progress has been made in the last two years in the measurement and QCD analysis 
of spin structure functions. The second generation CERN and SLAC experiments have been 
augmented by the HERMES experiment at DESY and these experiments now provide detailed 
information on the proton, neutron and deuteron spin structure. Here, we discuss measurements 
of the longitudinal asymmetry 



All 



where a^^ (a^^) is the cross-section when the lepton and nucleon spins are parallel (antiparallel) 
to each other. This asymmetry takes into account the beam and target polarisations as well as 
the dilution factor, the fraction of polarisable nucleons in the target. The asymmetry can then 
be related to the virtual-photon nucleon asymmetries Ai and A2 

All = D{A,+7]A2) 

where D is the depolarisation factor depending upon R = cl/o't and rj = t]{x, Q'^) is a kinematic 
factor. The spin structure function 

where 7 is a further kinematic factor and gi is then calculated using the known F2 and R values. 
The Ai term dominates and to a reasonable approximation the extracted structure function is 
proportional to the measured asymmetry 

^ 2x{l + R) D 



In the Quark Parton Model 



^?i(^,Q') = ^Ee?-AS(x,Q^) 



where AS = Y^tiiol ^ 0.1) + {ol ~ 0.1)] is the singlet summed quark and anti-quark distributions 
where {q^) is the quark distribution with the spins parallel (antiparallel) to the nucleon spin. 
In the QPM gi is therefore the charge-weighted vector sum of the quark polarisations in the 
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nucleoli. Similarly the non-singlet quark distribution Aq^s = J2i '^'^^t^^ [{qj — q^) + {qj — gf)], 
where <e^>= 2/9 for 3 light flavours is deflned such that this contribution evolves separately 
from the gluon {Ag) and singlet (AE) contributions when QCD corrections are taken into 
account. The DGLAP formalism outlined earlier can therefore be applied and Aq^s constrained 
by gi via QCD corrections. 

The corresponding integrals Fi = Jq gi{x, Q^)dx determine the total spin carried by the quarks 
at the measured Q^. In the QPM, below charm threshold 

14 1 1 
r1 = -i-Au + -Ad + -As) 



and, by isospin symmetry 



2'9 9 9 



11 4 1 
r" = -(-Au + -Ad+ -As) 
^2^9 9 9 ^ 



Beyond the QPM, the total spin of the nucleon can be written as the sum of the contributions 
from its constituents 

^ = i AE + + + Lg 

where the A's correspond to the intrinsic spins and the L's correspond to the angular momentum 
of the quarks and gluons. The proton spin puzzle is that only a fraction of the total spin is 
due to quarks. This puzzle remains, but flrst steps have now been taken to constrain the gluon 
contribution via the scaling violations of gi{x) and the valence and sea quark contributions via 
semi-inclusive measurements. 

The world data on gi{x) for the proton, deuteron and neutron are summarised in Fig. |11] 
for the CERN, SLAC and DESY experiments. [0 Comparison of the lower and higher 



data indicates that the scaling violations are relatively small over the well-measured range of 
1 < Q^^IO GeV^. The data are taken using a variety of targets (e.g. H(p), D{d) and ^He(n) 
for HERMES; NH3(p), ND3(rf), ^UD{d) and ^}ie{n) for the SLAC experiments; and, NH3(p), 
butanol(p) and deuterated-butanol((i) for SMC) and widely varying experimental techniques. 

There is no evidence for a rise of gf at the smallest values of x ~ 10-2 explored by SMC and 
E155. Indeed the QCD flt prediction is that gi becomes negative for x^lO~^ due to the relatively 
large positively polarised gluon at higher x driving the polarised gluon negative at small x. 
Clearly this will be of interest in the light of the HERA unpolarised results and of signiflcance 
in the extrapolations required for the determination of the sum rules. SMC presented data on 
the virtual-photon proton asymmetry, at low <Q^>= 0.01 GeV^ which extend to very low 
X = 10~^, as shown in Fig. |12|. [|35| The low x data were obtained with a dedicated low-Q^ 



trigger, requiring an observed hadron in each event which rejects radiative and other events 
with low depolarisation factors. At these very low values the data indicate that extreme 
QCD behaviour of gi oc l/xln^x (full line) proposed in |]5B[ is ruled out. However the less 
extreme QCD behaviours t^j' oc Inx or 2 + Inx indicated by the dotted lines give reasonable 
descriptions of the data. Clearly higher data are desirable in order to test the pQCD models, 
but the data do constrain the Regge behaviour. 

The values of gi in Fig. |TT] fall monotonically with decreasing x. The precision of this data 
is now approaching that of the proton data, which is important in the context of the Bjorken 
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Figure 11: Published world data on gi{x) from CERN experiments (0.003 < x < 0.6, 1 < 
< 100 GeV^) in the left column and SLAG and DESY experiments (0.014 < a; < 0.7, 
1 < < 20 GeV^) in the right column with statistical errors only. Results from two (MS 
scheme) NLO QGD fitting programs discussed in the text are superimposed. 
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Figure 12: SMC virtual-photon proton asymmetry, A^, for < >= 0.01 GeV^ data as a 
function of x compared to the behaviours discussed in the text. The systematic errors for the 
low {x, Q"^) and the reference set of data are indicated by the shaded bands. 

sum rule discussed below. All data are observed to be in good agreement, where the systematic 
errors (not shown) are typically smaller than the statistical errors. In addition, E155 presented 
new preliminary data on g\ and (?" (not shown) at = 5 GeV^ which are in good agreement 
with these pubhshed datasets. These data have very small statistical errors and extend the 
X range compared to the E143 data. This improved precision also requires that nuclear effects 
due to the assumed superposition of D and '^He states in the ^Li target data are understood, 
but possible uncertainties can be tested by comparison with data from different targets and 
other experiments. 



SMC have performed a NLO QCD fit to extract the singlet {x ■ AS) and gluon (x ■ /S.g) as 
well as the proton {x ■ Aq^g) and neutron (x ■ Aq^g) non-singlet polarised parton distributions 
shown in Fig. ^ Here, the deuteron spin structure function gf is assumed to be related to the 
proton and neutron structure functions by 

P , n _ ^^1 

where uj£, = 0.05 ± 0.01 is the D-wave state probability in the deuteron. 

The results of the fits from two NLO QCD programs are shown by the full and dashed lines in 
Fig. ITT| at the measured of each of the datasets. The comparison indicates that a good fit is 



obtained to the world data with x^/DoF = 127.4/(133 — 8), considering statistical errors only. 
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Figure 13: SMC QCD analysis of world gi data: singlet, gluon, proton and neutron non-singlet 
polarised parton distributions in AB scheme at = 1 GeV^. The fits with the statistical un- 
certainty are indicated by the cross hatched upper bands. The (small) experimental systematic 
uncertainties are indicated by the central bands and the (larger) theoretical uncertainties by 
the lower bands. 



The parton densities are parameterised and determined in analogy to the unpolarised case via 
the NLO DGLAP splitting kernels. Results are quoted in the Adler-Bardeen (AB) scheme, a 
modified version of the more conventional MS scheme, defined such that AEab is independent 
of Q^. These renormalisation/factorisation schemes are related via 



AE 



AB 



nf— A^ab(Q ) 

ZTT 



where ao(Q^) is the singlet axial-current matrix element and the A's correspond to the parton 
densities integrated over x. This scheme dependence is large and sufficient to create a negative 
AS(x) at small x and a smaller Ag{x) in MS scheme. However, the physical structure 
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functions and their integrals, such as ao(Q^) are unaffected by this definition. In Fig. 13, the 
singlet contribution is seen to be well constrained by gi. This remains true when account is 
taken of the theoretical uncertainties which include the variation of the renormalisation and 
factorisation scale by factors of two, varying as = 0.118 ± 0.003 within the given limits and 
changes to the starting scale and the functional forms of the parameters. Similarly the non- 
singlet contributions are reasonably well constrained and the rise of Ag^^ is mirrored by the 
fall of AgJ^g. However, Ag is rather poorly constrained especially when taking into account the 
theoretical uncertainties. Integrating over a; at = 1 GeV^ 

Ag = 0.99lJiI(stot.)l°:i(s2/s.)lJ:i(tWi/) 

in AB scheme compared to Ag = 0.25to^2i^tat.) in MS scheme. Consistent results are obtained 
in either scheme for the singlet contribution expressed in terms of the axial-current matrix 
element 

aoiQl = 1 GeV^) = 0.23 ± 0.07{stat.) ± 0.19{sys.) 

which compares to the QPM expectation ^ 0.58 and corresponds to about one third of the 
nucleon spin being carried by quarks. 

In the latest analysis by J. Ellis and M. Karliner, the world average value of the singlet 
matrix element in MS scheme is given by 

AS = 0.27 ±0.05. 

This approach utilises the 0{a^) calculations and 0{ag) estimates discussed below in relation 
to the Bjorken sum rule. Here, the consistency of the data taken at different values improves 
as successive higher-order QCD corrections are taken into account. 

Bjorken Sum Rule: The Bjorken sum rule is a fundamental prediction of QCD determined by 
the difference of the spins carried by the u and the d quarks 



r? - = / {gf - g-)dx = -f-\- C^iQ 

Jo O Qv 



6 gv 

where gA and gv are the axial-vector and vector weak coupling constants of neutron /?-decay 
and Ci^{Q'^) is the non-singlet perturbative QCD correction which has been calculated up to 
0{a^g) and estimated up to Oi^a^) using the same approach as for the GLS sum rule, [l^, |ll[ 
Here, higher twist contributions to the difference of the structure functions are assumed to be 
negligibly small. 

In the SMC NLO fits discussed above, the sum rule is imposed as an input with \ ^\ fixed rather 
precisely from /3-decay experiments. However, the NLO fits provide two methods to check the 
Bjorken sum rule. First, the input can be relaxed and ||^| fitted as an additional parameter in 
the global fit yielding 

r? - = 0.174 ± Qmb{stat)t'Zl{sys.)t',^,{theoTy) 

at g2 = 5 GeV^ which is in agreement with the theoretically expected value of 0.181 ± 
0.003. Second, the non-singlet contributions g\ — gi can be evolved to a common in order 
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Figure 14: Spin sum rules Ti versus F". Proton (EMS and SMC) and neutron (E142, E154 
and HERMES) data are indicated by the vertical and horizontal bands, respectively. Deuteron 
(E143 and SMC) data are indicated by the falling diagonal bands. The theoretical expectation 
from the Bjorken sum rule (= Ff — F") is indicated by the rising diagonal band while the 
Ellis- Jaffe sum rule is the area on this band indicated by the arrow. 
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Figure 15: Spin integrals, Fi proton, deuteron and neutron data versus compared to the 
Ellis- Jaffa sum rule (As = and SU(3)/ symmetry) expectation indicated by the shaded band. 
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to determine the sum rule more directly and with a minimum number of parameters. This 
yields a consistent value of 0.181'^Q^li{total) when evaluated at 0{a'^). This second method is 
potentially more precise, but awaits the high statistics data on Qi from E155 combined with 
the existing data on g"^ from E154. 

Ellis-Jaffe Sum Rule: Assuming that strange quarks do not contribute to the nucleon spin and 
SU(3) f symmetry, Ellis and Jaffe derived independent sum rules for the proton and neutron 

r? = -(-Au + -Arf) -0.17 

= -(-Au + -Ad) ~ -0.02 
^2^9 9 ^ 



in the high-Q limit and modified by singlet and non-singlet QCD corrections. In Fig. IT^ the 



world data on the Bjorken and Ellis-Jaffe sum rules are depicted graphically. As noted in 
relation to the SMC analysis, the data are in agreement and consistent with the Bjorken sum 
rule with a precision of around 10%. The origin of the spin puzzle was the EMC measurement 



of In Fig. |T5|, the world data on F^^, F^ and F" all indicate that the Ellis-Jaffe sum rule is 



broken at the 2-3a level. The strange sea quarks and/or the gluon therefore carry a significant 
fraction of the spin. It is currently impossible to distinguish an SU(3)/ symmetric sea or a 



As = (large gluon) solution in the NLO QCD fits. |38| A natural assumption would be that 
SU(3)/ symmetry is violated at the same level as in the unpolarised structure functions and 
the rest of the spin can be attributed to a large gluon polarisation, but this requires further 
experimental input. 



Semi- inclusive Asymmetries: SMC [39| and HERMES have recently produced data tagging 
the charge of final state hadrons. The asymmetry is defined in analogy to the inclu- 

sive case where the distribution of charged hadrons with momentum fraction z is statistically 
correlated with the struck quark flavour in a flavour tagging analysis. In Fig. |16| the extracted 
spin contributions Am„, Ady and Aq (introduced for x < 0.3) from SMC and HERMES are ob- 
served to be in agreement. Slightly different assumptions are made with respect to the strange 
quark sea where SMC assume an SU(3)/ symmetric sea and HERMES assume this symmetry 
is violated at the same level as in the unpolarised structure functions. However these effects are 
negligible for the predominantly pion final states in this LO QCD analysis. Both experiments 
observe Auy is positive and Ady is negative. The SMC results integrated over x are 

Auy = +0.77 ± 0.10(stat.) ± 0.08{sys.) 
Ady = -0.52 ± O.Uistat.) ± Om{sys.) 
Aq = +0.01 ± OM{stat.) ± 0.03{sys.) 

in agreement with the expectations from the NLO QCD fits. The polarised sea is compatible 
with zero although there are indications from the HERMES data that this contributes positively. 



Outlook: The analysis of polarised structure function data enables first measurements of the 
scaling violations to be performed. These provide a test of QCD and indicate that the gluon 
polarisation is positive. Similarly, the semi-inclusive measurements give first constraints on 
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Figure 16: Data from SMC and HERMES for the valence {Au^ and Ad^) and sea (Aq) spin 
densities versus x extracted from semi-inclusive asymmetry measurements. The dominant 
statistical errors only are shown. 
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the sea. The structure function data from E155 and semi-inclusive data from HERMES will 
soon provide further input. These data are, however, insufficient to determine the various spin 
parton distributions within the nucleon. By analogy with Table ^ we can therefore consider a 
future programme of spin structure measurements which will enable the partonic spin structure 
to be unravelled, as shown in Table 0. In the next few years HERMES, COMPASS and the 
RHIC experiments will focus on the determination of the gluon. It is also possible that the 
polarised proton technology developed at RHIC could be utilised at HERA in order to explore 
polarised structure functions at low {x,Q^). With a high luminosity HERA machine the data 
would also extend to high where neutral current and charged current events would provide 
detailed information on this spin structure. 



4 HERA Results 

The HERA ep collider has improved its performance in successive years, providing large e~^p 
datasets from 1994-97 running. During this period 27.5 GeV positrons collided with 820 GeV 
protons resulting in data samples of 37pb~^ and 46pb~^ for HI and ZEUS, respectively. As 
a measure of the progress which has been made in the study of the scaling violations of F2, 



Fig. [T^ illustrates how the latest HERA data, presented for the first time at ICHEP98, extend 
the reach in beyond 10,000 GeV^ with x extending up to 0.65. At high values, the effects 
of Z° exchange are significant and the value of F2"^ is quoted where 

p _ pem I Q pint , ( Q_ ] ™fe 

where -^2^"^, -F2"* contributions due to photon exchange, interference and 

Z exchange, respectively. The data also extend to very low x below 10~^: here it should be 
noted that the precision of the low-a;£0.03 data is now comparable to that of the fixed-target 
data at higher x. There is a region of overlap where the HERA and fixed-target experiments 
can be compared. In particular, the analysis of the ZEUS results at low y may resolve the 
CCFR-NMC discrepancy discussed earlier in the context of the "5/18ths rule". However, this 
comparison (not shown) is currently inconclusive, with the ZEUS data lying between the NMC 
and CCFR data. 



4.1 Low-Q^ Results 

Transition Region: The rise of F2 with decreasing x or, equivalently, the rise of o"**^ with 
increasing W has stimulated significant experimental and theoretical developments in the un- 
derstanding of QCD at high energies. One challenge is to explore how and where the transition 
occurs from soft to hard physics and interpret low-Q^ data. Measurements have been performed 
using dedicated low-angle taggers (e.g. the ZEUS BPC) and shifted vertex (SVX) data |^ as 
well as QED Initial State Radiation (ISR) data |0 in order to map out this region. In Fig. |l^. 



a compilation of the latest measurements available from HERA and E665 indicates that the 
different experiments and techniques agree with a precision of around 5% from the most recent 
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data. The significant rise of F2 is apparent for Q'^>\ GeV^, a behaviour which is described by 
the ZEUS NLO (pQCD) fit. This behaviour is not reproduced in the DL (Donnachie-Landshoff 
Regge) fit. The ZEUS collaboration has performed both types of fits to the F2 data exploring 
this transition region. ||4l| In order to determine the relationship between low-Q^ ZEUS BPC 



data measured in the range 0.1 < < 0.65 GeV^ and = GeV^ data, a Generahsed 
Vector Meson Dominance (GVMD) approach can be taken. GVMD relates the virtual-photon 
cross-section to the real cross-section via 

^tot _ tot ^^^0 

for fixed W and (Tl = 0. A good description of the data is found with Mq = 0.53 ± 0.04 ± 
0.10 GeV^. Regge theory then determines the W dependence of the data, combined with lower 
energy photoproduction experiments, as 

where the Reggeon intercept is fixed to 0.5 from hadropro duct ion data and lower energy 
photoproduction data also constrain Ap, A^i and a p. From the BPC data alone, the pomeron 
intercept value is 

^BPC ^ ^_^4^ ^ 0.020(stot.) ± O.OU{sys.) 

to be compared with the Donnachie-Landshoff value = 1.08. In this range, the rise of 
the cross-section is therefore relatively modest. Combining the GVMD and Regge approaches, 
the resulting ZEUSREGGE fit is used to parameterise the and W dependence of the low 
g2 < 1 GeV2 data. The ZEUS NLO QCD fit to the ZEUS 94 and ZEUS SVX data, incorpo- 
rating NMC and BCDMS data, (discussed below) is used to determine the behaviour and the 
uncertainties on the gluon and singlet quark densities at low x for > 1 GeV^. 

To quantify the behaviour of F2, fits to the E665 and ZEUS data of the form F2 = c - x'''^"^ \q2 
are performed. The parameter X^s — Q^p — 1 for x < 0.01 is then plotted as a function of in 
Fig. [TP|. A relatively slow transition from Aofr 0.1 is observed with increasing Q^. Also shown 
are fits to the DL and ZEUSREGGE parameterisations, fitted over the same x range as the 
data, for < 1 GeV^. These describe the data reasonably but are systematically lower. For 
> 1 GeV^, the data are compared to the GRV94 prediction where the starting scale for the 
evolution of the parton densities is rather low (~ 0.3 GeV^) and pQCD evolution generates the 
rise at small x: this approach is observed to reasonably describe the data and the description 
is improved using the GRV98 PDF (not shown). The rise at small x is also described by the 
ZEUSQCD fit, where the ZEUS data is used as an input. 

This rise of F2 with decreasing x is intimately coupled to the scaling violations via the gluon 
density (in leading order dF2/d\iaQ'^ ~ xg{x) neglecting sea quark contributions). In Fig. 
fits of the form F2 = a + b ■ ln{Q'^)\ri. have been performed to the ZEUS data and the param- 
eterisations discussed above. For x ~ 2.10"^, corresponding to < >~ 4 GeV^ there 
is a qualitative change in behaviour where the scaling violations stabilise and then decrease 
for lower-x values, a behaviour which is not reproduced by the GRV94 PDF. The question is 
whether this scaling violation behaviour and the slow onset of the rise of F2 with decreasing x 
can be simultaneously understood. 
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A. Mueller has discussed the scaling violation behaviour in terms of a geometric model where 
the spatial extent, Rq, of the qq fluctuation of the virtual photon in the 7*p fixed-target frame 
is related to the height of the plateau in the scaling violations p3| 



dF,, 



rflnQ2 



Ana 



For 



dF2 



0.4 GeV ^ and ^ 4 GeV^, the spatial extent Rq ^ 0.3 fm. This appears to 
e at which a transition takes place and the partons in the proton start to overlap. |3] 



dln( 

be the sea . 

However, perturbative QCD is pushed to its limit and it will be important to test that the 
parton densities extracted from F2 can be universally applied. 



The ZEUS NLO fit to the > 1 GeV^ data describes the data, demonstrating that there is 
sufficient flexibility in such an approach to go down to relatively low Q^. However, the relatively 
stable scaling violations observed around <Q^ >~ 4 GeV^ in Fig. yield a gluon contribution 
which is rapidly diminishing at small-x and which is in fact smaller than the singlet sea quark 
contribution for small starting scales. For larger values the gluon dominates the sea and we 
have an intuitively appealing picture where gluons radiate sea quarks whereas, in this low-Q^ 
region, the sea appears to be driving the gluon at low x. Whether such low-Q^ partons are 
universally valid could be tested using e.g. low-Q^ diffractive vector meson data. R H3 



An important part of the ZEUS NLO QCD fit is the determination of the uncertainties on 



the gluon and singlet quark densities. These are given for the gluon distribution in Fig. 
The overall uncertainty is estimated by combining in quadrature: the experimental systematic 
uncertainties on the ZEUS as well as the NMC and BCDMS data; the theoretical uncertainties 
on as(M|) by ibO.005, the relative strange quark content by ±50% and the charm mass by 
±0.2 GeV; and, the parameterisation uncertainties on the starting scale, by varying 1 < < 
7 GeV^, as well as using a more flexible form of the input in terms of Chebycheff polynomials 
and redefining the including stat. © sys. errors. These variations correspond to a precision 
on the gluon of ~ 10% at = 20 GeV^ where the renormalisation and factorisation scales 



are set to ii\ = = Q"^ . The role of the scale uncertainties is discussed in with respect 
to future capabilities to determine as(M|) at HERA. The theoretical and parameterisation 
uncertainties are amplified at low such that the gluon is rather poorly determined from the 
scaling violations of F2 in the transition region discussed above.^ It is clear, however, that the 
gluon is significantly suppressed at low Q"^ . 

F2 Determination: D* {Ktt)7[s measurements in DIS provide a significant test of the gluon 
density of the proton determined from the scaling violations of F2. They also help to constrain 
theoretical uncertainties in the fits to F2 where different prescriptions for charm production are 
adopted. The ZEUS preliminary cross-section ct^p^-D'x _ 3.55 j^o.31]'^o;5o ^t) is measured in the 
range 1 < < 600 GeV^, 0.02 <y < 0.7, 1.5 < p^* < 15 GeV, and '{rj^'l < 1.5. In general, 
the HI 1^ and ZEUS data agree with the Harris-Smith NLO calculations where the 



fraction /(c D*+) = 0.222 ± 0.014 ± 0.014 is determined from LEP data [||], the Peterson 
fragmentation function is characterised by ec = 0.035 and the renormalisation and factorisation 



''^ The gluon can even become negative at low (in MS scheme) but the physical quantities, F2, F| and 
Fl remain positive within the quoted errors. 
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scales are set to /x^ = = Q^ + 4m^. There is, however, a small discrepancy at lower xr,* and 
higher 77^* , corresponding to the proton direction, where the data lie above the prediction. A 
similar discrepancy is also observed in the first analysis in the semi-leptonic decay channel. 



Together these results indicate that the fragmentation of charm in ep processes is worthy of 
further investigation. At this stage, it is reasonable to extrapolate the measured cross-section 
to the full {ti^Pt} range0 to determine -^2(0;, Q"^) via the expression 



dx dQ"^ xQ^' 



■[Y^F^{x,Q')-y'Fl{x,Q')] 



where the contribution can be estimated as a QCD correction. In Fig. ^ the HERA 
F2{x, Q^) data mirror the rise of F2 at small x. The data are in agreement with the GRV94 
PDF, where the band represents an estimated theoretical uncertainty due to the effective charm 
mass [rric = 1.4 ± 0.2 GeV). This comparison verifies the steep rise of the gluon density at low 
X with a precision of ^ 15 — 20%. 

Fl Determination: The contribution of F^ enters as a QCD correction to the total DIS cross- 
section where F^ = F2 — 2xFi. As such it provides an additional method to calibrate the 



gluon at low x. HI |42] have used two methods to extract F^ from the reduced cross-section 



2 

a = F2 — ^ ■ Fi a.t high y. This is the region where the scattered electron energy is low: in the 
HI analysis scattered positrons are measured down to E^' > 3 GeV and backgrounds reduced 
by requiring the associated track to have correct charge. F^ is determined as a function of 
> 3 GeV^ by measuring local derivatives of da /d logy and observing deviations from a 



straight line at high y. These data are denoted by the stars in Fig. [23. Here, the data 



are compared to an earlier extrapolation method |^ applied to the same data, which yields 



consistent results, as well as to the HI NLO QCD fit to HI, NMC and BCDMS F2 data. The 
data are in agreement with the QCD expectations although there is an indication of a relatively 
large Fl contribution at the highest y (corresponding to lowest x) values. In conclusion, a 
consistent value for the gluon density at low x£10~^ may be extracted from the data on F2, -F| 
and Fl with a precision of ~ 10% at = 20 GeV^. 



4.2 High-Q^ Cross-Sections 

The HERA collider provides a unique window to explore ep interactions at the highest energies, 
extending the range of momentum transfer by about two orders of magnitude compared to 
fixed-target experiments. As the HERA luminosity increases we explore the region of ~ 
10^ GeV^, where electroweak effects play a role. It is in this unexplored kinematic region that 
we are sensitive to deviations from the standard model (SM). 



In 1997, HI and ZEUS reported an excess of events compared to the SM predictions 
from the neutral current (NC) data taken during 1994 to 1996. For the HI analysis an accu- 
mulation of 7 events in a reconstructed e+g mass window of 200±12 GeV was found, compared 
to an expectation of 0.95 ± 0.18 from 15 pb~^ of data. One further event was found from the 



^^This procedure neglects the possibility of additional contributions outside the measured region due, for 
example, to intrinsic charm. 
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1997 data corresponding to a further 22 pb~^, yielding 8 events compared to an expectation of 
3.01 ± 0.54. For the ZEUS analysis the observed rates agreed with expectations except for an 
excess at the highest where two outstanding events with ~ 40, 000 GeV^ were observed 
from a luminosity of 20 pb~^. These events still clearly stand out but no new NC outstanding 
events are observed in the 1997 data, corresponding to a further 26.5 pb~^. Similarly in the 
charged current (CC) channel, the number of events is higher than expectations but is consis- 
tent with the standard model. Attention has therefore focussed on measuring the cross-sections 
at the highest accessible Q"^ values. 

The theoretical uncertainty on the cross-sections was determined as discussed w.r.t. the ZEUS 
low-Q^ NLO QCD fit using high-x Ff" and F2 data to yield SM cross-section uncertainties of 
~ 6-8% on the NC cross-section and ~ 6-12% on the CC cross-section at the highest accessible 
values. These cross-sections therefore represent a benchmark for the standard model. The 
cross-sections, discussed below, are corrected to the electroweak Born level and integrated over 



the measured range of y for HI and corrected to the complete y range for ZEUS. 



Neutral Current Cross-Sections: High-Q^ neutral current events are easily identified from the 
high-energy scattered positron. The cross-section is particularly sensitive to the valence M-quark 
distribution in the proton 

Here, F2 is the generalised structure function, incorporating 7 and Z terms, which is sensitive 
to the singlet sum of the quark distributions {xq + xq) and xF^ is the parity-violating [Z- 
contribution) term which is sensitive to the non-singlet difference of the quark distributions 
[xq — xq). The data are now becoming sensitive to electroweak 7Z interference effects, which 
suppress the NC cross-section by ~ 30% for > 10,000 GeV^. 



In the upper plot of Fig. the HI cross-section is observed to fall over more than six orders 
of magnitude. The ratio of the data to the SM, adopting the HI NLO QCD fit, is shown in 
the lower plot of Fig. p4| where agreement is observed up to ~ 30,000 GeV^. Comparison 
of the the data uncertainties with those from theory (not shown) indicates that the data will 
constrain the parton densities of the proton at large x. 

Charged Current Cross-Sections: Charged current events are identified by their missing trans- 
verse momentum (p^) due to the escaping neutrino. The cross-section is sensitive to the valence 
d-quark distribution in the proton 



dV,+„ Gl ( M2 



w 



[u + c+ {l-yf{d + s)]. 



In the ZEUS analysis, da'"'" /dQ"^ was measured for > 400 GeV^ using the Jacquet-Blondel 
method where Qj^ = p'^/l—yjB, with an RMS resolution on ~ 25%, reflecting the 35%/ Ve 
hadronic energy resolution. The systematic uncertainties, mainly due to the hadronic energy 
scale uncertainty of ±3%, correspond to ~ 15% uncertainties on the cross-section at lower Q'^ 
but increase at larger x and Q^. 



29 



In the upper plot of Fig. ^ the cross-section is observed to fall over more than four orders of 
magnitude. The ratio of the data to the SM, adopting the CTEQ4D PDF, is shown in the lower 
plot of Fig. where good agreement is observed up to of ~ 10,000 GeV^. Comparison of 



the the data uncertainties with those from theory (shaded band) indicates that the data will 
help to constrain the d-quark densities at large-x. The cross-section is suppressed relative to 
NC 7 exchange due to the propagator term: this characteristic dependence on has been 
fitted to yield a value for the mass of the exchanged (space-like) W^-boson of 

Mm/ = 78.6l|^(stot.)+;^:|^(si/s.) GeV (ZEUS prelim.) 
Mw = 81.2 ± 3.3{stat.) ± 4.3(s?/s.) GeV (HI prelim.) 

with an additional PDF uncertainty of ±1.5 GeV estimated in the ZEUS analysis and an 
uncertainty of ±3 GeV due to electroweak radiative corrections estimated in the HI analysis. 

Returning to the Bodek-Yang analysis, fl^ discussed in relation to the fixed-target results, the 
ZEUS CC cross-section is plotted as function of x for > 400 GeV^ in Fig. |26|. The increase in 



the ratio of d/u corresponds to an increase of the CC cross-section at high x. The uncertainties 
on the data are large in this region, but this modification does result in better agreement with 
the data than the standard PDF's at large x. 

A comparison of the HI and ZEUS NC and CC cross-sections for > 1000 GeV^ is given 
in Fig. ^ The NC and CC data are in agreement and both cross-sections agree with the 
SM prediction over a broad range range of Q^. At high > 10,000 GeV^, the CC cross- 
section is suppressed relative to the NC cross-section due to the d/u ratio being less than unity. 
The measurement of the HERA CC/NC ratio of cross-sections will therefore provide a direct 
determination of this ratio, free from the uncertainties due to nuclear binding effects. 



5 Summary and Outlook 

There were many highlights in our deepening understanding of nucleon structure presented 
at the ICHEP98 conference. CCFR have performed final analyses of their structure function 
data which lead to precise tests of QCD. NuTeV already provide electroweak input and aim 
to reduce the uncertainty on as{M^) to ±0.002. Important developments have been made in 
the understanding of higher twist effects in terms of renormalon theory. Input from various 
experiments and a reassessment of the importance of nuclear binding effects in the deuteron 
lead to the conclusion that the ratio of d/u parton densities increases, whilst the ratio of d/u 
decreases compared to earlier determinations. 

The second phase of spin experiments are now exploring spin structure via NLO QCD fits. 
The first observations of scaling violations from combined fits to the world data indicate that 
the source of the spin puzzle hes with the gluon spin density, Ag, although large uncertainties 
remain. Semi-inclusive measurements now provide input on the quark composition of the spin. 

The first discovery at HERA was the rise of F2 at low x. Precise data now enable the rise of 
the associated scaling violations with decreasing x and hence the gluon to be determined. The 
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determinations of F2 and enable this gluon distribution to be calibrated. The fall of these 
scaling violations at low (x, Q"^) enables the region where parton confinement effects take place 
to be explored at high energies. 

The large e^p data sample enables cross-sections to be measured at very large where elec- 
troweak effects start to play a role. The HERA data are consistent with the Standard Model 
and place constraints on the parton densities at large x. The outlook is for a similar sample 
of ~ 50pb~^ of e~p data in the next two years, prior to the HERA upgrade where luminosities 
will be increased five-fold. Deep inelastic scattering has historically led to the discovery of the 
nucleus, quarks and electroweak neutral currents. The discovery potential of current and future 
experiments is high and the field continues to provide important input to our understanding of 
sub-nuclear structure. 



Acknowledgements 



It is a pleasure to thank Alan Astbury and co-organisers for an excellent conference. Many 
thanks to Halina Abramowicz, Arie Bodek, Antje BruU, Allen Caldwell, Abe Deshpande, John 
Dainton, Robin Devenish, Martin Erdmann, Laurent Favart, Thomas Gehrmann, Tim Green- 
shaw, Beate Heinemann, Peppe lacobucci, Robert Klanner, Max Klein, Masahiro Kuze, Ludger 
Lindemann, Alan Martin, Gavin McCance, Jan Okrasihski, Jen-Chieh Peng, Alex Prinias, 
Robert Waugh, Arnulf Quadt, David Saxon, Stefan Schlenstedt, Mike Vetterli, Manuella Vinc- 
ter, Bryan Webber, Jim Whitmore, Un-Ki Yang, Rik Yoshida, Jaehoon Yu and all the speakers 
in the structure functions parallel session for their help, inspiration and advice. I am grateful 
to the Alexander von Humboldt Foundation, DESY and PPARC for their financial support. 



References 

[1] V. Fadin and L. Lipatov, |hep-ph/9802290 . 



[2] J.R. Forshaw, G.P. Salam and R.S. Thorne, |hep-ph/9812304 



[3] L.V. Gribov, E.M. Levin and M.G. Ryskin, Phys. Rep. 100, 1 (1983); E.M. Levin and 
M.G. Ryskin, Nucl. Phys. B 18C, 92 (1990). 

[4] J. Huston, ICHEP98 proceedings. 

[5] Y. Dokshitzer, ICHEP98 proceedings. 

[6] M. Erdmann, ICHEP98 proceedings. 

[7] A.D. Martin, R.G. Roberts, W.J. Stirling and R.S. Thorne, |hep-ph/9803"445 
M. Gliick, E. Reya and A. Vogt, |hep-p/980640i . 



[9] CTEQ Collab., H.L. Lai et al., Phys. Rev. D 55, 1280 (1997). 

31 



[10] S.A. Larin and J.A.M. Vermarseren, Phys. Lett. B 259, 345 (1991). 

[11] A.L. Kataev and V.V. Starshenko, Mod. Phys. Lett. A 10, 235 (1991). 

[12] J.H. Kim, D.A. Harris et al., |hep-ex/98080T5| , submitted to Phys. Rev. Lett. 



[13] S.A. Kulagin, |hep-ph/ 9809219 



[14] A.L. Kataev, G. Parente and A.V. Sidorov, |hep-ph/ 9809500 



[15] CCFR Collab., W. Seligman et al., |hep-ph/97010T7 



[16] M. Virchaux and A.M. Milsztajn, Phys. Lett. B 274, 221 (1992). 

[17] C. Caso et al.. Euro. Phys. C 3, 1 (1998). 

[18] F. Eisele, Physics in Collision proceedings, |liep-ph/9807028 . 



[19] U. Yang and A. Bodek, ICHEP98 proceedings, |liep-ph/9809480 



[20] L. Frankfurt and M. Strikman, Phys. Rep. 160, 235 (1988). 
[21] M. Dasgupta and B.R. Webber, Phys. Lett. B 382, 273 (1996). 



[22] CCFR Collab., U.K.Yang et al., DIS98 proceedings, |hep-ex/ 980602^ ; J. Yu, ICHEP98 
proceedings. 



[23] CDF Collaboration, F. Abe et al., |hep-ex/9809"001| , submitted to Phys. Rev. Lett.; W. 
Sakumoto, ICHEP98 proceedings. 

[24] NMC Collab., M. Arneodo et al., Phys. Rev. D 50, 1 (1994). 

[25] S. Kumano and K. Umekawa, |hep-ph/9803359 . 



[26] E866 Collab., E.A. Hawker et al, Phys. Rev. Lett. 80, 3715 (1998). 

[27] E866 Collab., J.C. Peng et al., |hep-ex 980428^ , submitted to Phys. Rev. D. 



[28] HERMES Collab., K. Ackerstaff et al., |liep-ex/9807013| ; M. Vincter, ICHEP98 proceedings. 



[29] C. Boros, J. Londergan and A. Thomas, |hep-ph/9804410| . C. Boros, ICHEP98 proceedings. 
[30] CCFR Collab., A. Bazarko et al., Z. Phys. C 65, 189 (1995). 
[31] D. Karlen, ICHEP98 proceedings. 

[32] SMC Collab., B. Adeva et al., CERN-EP/98-86, submitted to Phys. Rev. D; A. Deshpande, 
ICHEP98 proceedings. 

[33] G. AltareUi, R.D. Ball, S. Forte and G. Ridolfi, Nucl. Phys. B 496, 337 (1997). 

[34] L. Sorrell, ICHEP98 proceedings. 

32 



[35 
[36 

[37; 

[38 
[39 
[40 
[41 

[42 

[43 

[44; 

[45 
[46 

[47; 

[48 
[49; 
[50 
[51 

[52; 



E. Rondio, ICHEP98 proceedings. 

F. E. Close and R.G. Roberts, Phys. Lett. B 336, 1257 (1994). 

J. Ellis, DIS98 proceedings (World Scientific) in press; J. Ellis and M. Karliner, Phys. Lett. 
B 341, 397 (1995). 



B. Lampe and E. Reya, |hep-ph/9810"270 



SMC Collab., D. Adeva et al., |hep-ex/971100^ , submitted to Phys. Lett. B. 
C.A. Miller, ICHEP98 proceedings. 

ZEUS Collab., J. Breitweg et al., DESY 98-121, submitted to Euro. Phys. C; A. Quadt, 
ICHEP98 proceedings. 

HI Collab., C. Adloff et al., ICHEP98 contributed papers 534 and 535; M. Klein, ICHEP98 
proceedings. 

A. H. Mueller, DIS98 proceedings (World Scientific) in press. 
M. Botje, DIS97 proceedings, AIP Conf. Proa. 407, 393 (1997). 

HI Collab., C. Adloff et al., ICHEP98 contributed papers 538 and 540; E. Tzamariudaki, 
ICHEP98 proceedings. 

ZEUS Collab., J. Breitweg et al, ICHEP98 contributed papers 768 and 772; W. Verkerke, 
ICHEP98 proceedings. 

B. W. Harris and J. Smith, Nucl. Phys. B 452, 109 (1995); B.W. Harris and J. Smith, 
Phys. Lett. B 353, 535 (1995); B.W. Harris and J. Smith, |hep-ph/ 9706551 . 

OPAL Collab., K. Ackerstaff et al.. Euro. Phys. C 1, 439 (1998). 

HI Collab., S. Aid et al, Phys. Lett. B 393, 452 (1997). 

HI Collab., C. Adloff et al., Z. Phys. C 74, 191 (1997). 

ZEUS Collab., J. Breitweg et al., Z. Phys. C 74, 207 (1997). 

HI Collab., C. Adloff et al., ICHEP98 contributed paper 533; M. Fleischer, ICHEP98 
proceedings. 



[53] ZEUS Collab., J. Breitweg et al., ICHEP98 contributed papers 751 and 752; J. Grosse- 
Knetter, ICHEP98 proceedings. 



33 



Table 1: Processes studied in global parton distribution fits, -k indicates new preliminary data 
presented at ICHEP98, included in this review, f indicates final data presented at ICHEP98. 
* indicates data used in the MRST fits. (Courtesy of A. Martin.) 



Process/ 
Experiment 



Leading order 
subprocess 



Peirton behaviour probed 



DIS {fxN fxX) 

^2 ,^2 ,^2 1^2 

(SLAC, BCDMS, 
NMC, E665)* 

DIS {vN ixX) 

jpuN ^jpuN 
^2 )'^-'^3 

(CCFR)* 



t 



DIS (small x) 

F^^ (HI, ZEUS)* ★ 

DIS (Fl) 

NMC, HERA ★ 
£N ccX 

F| (EMC; HI, ZEUS)* ★ 

uN fi+iJL-X 

(CCFR)* 

pN jX 

(WA70*, UA6, E706, ...) 

pN n+n'X 

(E605, E772)* 

pp,pn — > ii^ix~X 

(E866, NA51)* f 

ep, en enX 
(HERMES) t 

pp WX(ZX) 

(UAl, UA2; CDF, DO) 

cisym (CDF)* f 

pp tiX 

(CDF, DO) 

pp ^jet + X 

(CDF, DO) 



W*q q' 
r{Z*)q ^ q 

1*9 qq 

7*c c 
W*s c 
qg 7g 

qq^i* 

uu, dd 7* 
ud, du — >• 7* 

7*g — > q with 
q — u,d,u,d 

ud^W 

qq, 99 tt 
99, q9, qq 2j 



Four structure functions 

u + u 
d + d 
u + d 

s (assumed = s), 
but only / xg{x, Ql)dx ~ 0.35 
and /(J— ■u)dx ~ 0.1 

A 

{xq ~ xg ~ x~^^) 

9 

c 

(a;> 0.01; a; < 0.01) 
s l{u + d) 

(y( at X ~ 2pj,/y/s — >• 
X ^ 0.2-0.6 



q — ...{1 — x) 



vs 



u-d (0.04 < a; < 0.3) 



u-d (0.04 < a; < 0.2) 



M, d at x ~ Mw I \/s — > 
x^0.13; 0.05 

slope oi u/d at X ^ 0.05 — 0.1 

q,g at X > 2mt/ \/s — 0.2 

g,fif at X ~ IE-tI^s 
X ^ 0.05 - 0.5 



34 



Table 2: Processes studied in global spin parton distribution fits. The first group of datasets 
correspond to current experiments. The second group correspond to near-future (year 2000) ex- 
periments. The third group correspond to potential year 2005 experiments at HERA. (Courtesy 
of T. Gehrmann.) 



Process/ 
Experiment 


Leading order 
subprocess 


Parton behaviour probed 


DIS (£N iX) 

^jid f~i^T^ 
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m h+h-X 

(COMPASS) 
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pp — > jets X 
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^9 (?) 
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(HERA) 


W*q q' 
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DIS (small x) 

gT (HERA) 


7*g q 


{Aq ~ x°''',Ag ~ x°'o) 


IN ^ £ jets X 
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i*q qg 
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(HERA-N) 
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Figure 17: Scaling violations of F2 versus for various x ranges from the latest HI and ZEUS 
preliminary data compared to fixed-target experiments. (The quoted x values correspond to 
the ZEUS measurements or the nearest x values from other experiments). 
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Figure 18: HERA F2 data for various intervals as a function of x exploring the transition 
region. The data are compared to the Donnachie-Landshoff Regge model (lower dotted line) 
and the ZEUS NLO QCD fit (full line) discussed in the text. 
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Figure 19: X^s versus and < a; > from fits to ZEUS and E665 data of the form F2 
c ■ x"'^"'^ \q2. The data are compared to the QCD and Regge fits discussed in the text. 
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Figure 20: dF2/d\nQ^ versus x and <Q^ > from fits to ZEUS data of the form F2 — a + b ■ 
ln((5^)|a;. The data are compared to the QCD and Regge fits discussed in the text. 
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Figure 21: Gluon momentum distribution, xg{x) at g2 = 1, 7 and 20 GeV^ in MS scheme 
including the uncertainties discussed in the text. The shaded band corresponds to the ex- 
perimental and theoretical uncertainties and the hatched band indicates the parameterisation 
errors discussed in the text. 
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Figure 22: HERA F2 data for various intervals as a function of x. The data are compared 
to the Harris-Smith NLO QCD calculation using the GRV94 PDF input discussed in the text. 
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Figure 23: HI determination of Fl versus for y = 0.68 and y = 0.82 using the derivative 
method (*) compared to the pubhshed method (•) and the NLO QCD fit expectation (shaded 
band). The upper hmit = F^ is indicated by the full line. 
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Figure 24: HI neutral current cross-section da'^^ /dQ"^ versus 

fo^ y < 0.9 and E^, > II GeV 

(upper plot) and ratio with respect to the standard model prediction (lower plot). The shaded 
band represents the luminosity uncertainty of ±2.6%. 
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Figure 25: ZEUS charged current cross-section da^"-^ /dQ"^ versus for < ?/ < 1 (upper plot) 
and ratio with respect to the standard model prediction (lower plot). The shaded band in the 
ratio plot represents the ZEUS NLO QCD fit uncertainty. 
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ZEUS CC Preliminary 1994-97 




Figure 26: ZEUS charged current cross-section da^'" jdx versus x for > 400 GeV^ and ratio 
with respect to the standard model prediction (lower plot). The upper curve in the ratio plot 
represents the Bodek-Yang d/u modified MRS(R2) parameterisation discussed in the text. 
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HERA e p DIS cross section 94 - 97 




Figure 27: HERA e~^p DIS cross-sections neutral and charged current data for y < 0.9 compared 
to the standard model prediction, adopting the CTEQ4D PDF. 
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